Pairbreaking Without Magnetic Impurities in Disordered Superconductors 
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We study analytically the effects of inhomogeneous pairing interactions in short coherence length 
superconductors, using a spatially varying Bogoliubov-deGennes model. Within the Born approx- 
imation, it reproduces all of the standard Abrikosov-Gor'kov pairbreaking and gaplessness effects, 
even in the absence of actual magnetic impurities. For pairing disorder on a single site, the T-matrix 
gives rise to bound states within the BCS gap. Our results are compared with recent scanning tun- 
neling microscopy measurements on Bi2Sr2CaCu208+i with Zn or Ni impurities. 
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The effect of disorder in superconductors has long been 
a subject of considerable interest. A generally accepted 
physical picture is that magnetic impurities destroy su- 
perconductivity by locally breaking the pairs, ||, ^, ^] 
whereas non-magnetic impurities are not pairbreaking, 
according to Anderson's theorem. |^ This is true for an 
isotropic s-wave BCS superconductor, in which the or- 
der parameter is uniform and momentum independent. 
Since in most high transition temperature Tc cuprates, 
the suppression of with Zn or Ni doping is comparable, 

there were proposals to explain this in terms of d-wave 
superconductivity. However, nuclear magnetic reso- 
nance (NMR) experiments indicated that the nominally 
non-magnetic Zn++ ions polarize the spin background 
in the Cu02 planes upon substitution of the nominally 
5=1/2 Cu++ sites. §, | Moreover, recent scanning tun- 
neling microscopy (STM) measurements directly above 
the Zn or Ni impurity sites observed strong resonance 
peaks. 1^, |l^ Very recently, several groups noticed from 
STM measurements that underdoped Bi2Sr2Cu208+5 
(BSCCO) and Bi2-3;Pba;Sr2Cu208-i-5 are extremely dis- 
ordered on a scale of a few nm. [ij, |l^ [m) This 
disorder is characterized by two gaps, one corresponding 
to the superconducting gap, with characteristic super- 
conducting peaks, and a non-superconducting gap. 

There is now a large body of evidence that the pseu- 
dogap observed in cuprates above Tc is not supercon- 
ducting. The most convincing of these experiments de- 
termined that the magnetic field dependence of the re- 
sistivity, NMR spin-lattice relaxation rate, magnetiza- 
tion, and the gaps observed in intrinsic tunneling ex- 
periments for the superconducting and pseudogaps are 
qualitatively different. ^ In particu- 
lar, the pseudogap regime is field independent until one 
reaches the Zeeman field for breaking up chargeless spin 
zero pairs. This is precisely consistent with the pseu- 
dogap being particle-hole pairs, which could be either 
of the charge-density wave (CDW) or spin-density wave 
(SDW) form. As in one dimension, such excitations are 
expected to arise from repulsive interactions between like 
charges, whereas superconductivity can arise when the 



interactions between like charges are attractive. 

Thus, if indeed the disorder involves a percolation 
problem between superconducting and density-wave re- 
gions on a nanometer scale, then the phase coherence 
can only arise from Josephson coupling of the supercon- 
ducting grains. In addition, one would expect the c-axis 
tunneling to be very incoherent, as inferred in BSCCO 
from c-axis twist Josephson junction experiments, pof 

In this letter, we assume the superconductor is elec- 
tronically disordered on the scale of the coherence length. 
We further assume the essential ingredient in the disor- 
der is not one of impurities, but rather disorder in the 
pairing interaction itself. Thus, we expect that Tc varies 
from site to site, as does the resulting order parame- 
ter amplitude. We treat this type of disorder using a 
Bogoliubov-deGennes procedure, assuming that the or- 
der parameter amplitude varies locally. j2^, 

Here we show that in the Born approximation, this 
problem maps exactly onto that of pairbreaking in a su- 
perconductor, with all of the features of that model, in- 
cluding gapless superconductivity. [Q, ||, |[ ^ At a partic- 
ular defect site, the T-matrix gives rise to bound states 
within the gap, even without magnetic impurities. 

We use the Nambu representation. 
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c|(r) 



Vcj(r)/ 

ft(r) = (c|(r), c|(r), c^{r), c^{r)) , (1) 

where C|(r) [c|(r)] annihilates [creates] a quasiparticle 
with spin eigenstate | [t] at the position r. We set h — 
c = ks = The Hamiltonian under study is H = 
Ho + Hi + H2 + H3, where 



Ho 



d^r¥{r)[i{r)p3(Jo + AoP2f72]*(r), (2) 



Vi(y) = Ui{r)p3ao, 



(3) 
(4) 
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Mr) = U2{r)S-l^/[S{S+l)]'/^ 
= xpafTi + ypocr2 + zpsCTs, 



(5) 
(6) 



where Pj<Tj' = Pj (E) CTj' is a rank 4 tensor composed of 
two Pauli matrices for j, j' — 1,2,3 and po: ctq are rank 2 
identity matrices, respectively. Hq is the BogoHubov- 
deGennes version of the BCS Hamiltonian, with mo- 
mentum space quasiparticle energy dispersion relative 
to the Fermi energy /i, Aq{T) is the real bare uniform 
BCS order parameter, and Hi and H2 are the inter- 
actions due to scattering off random non-magnetic and 
magnetic impurities with effective potentials Ui{r) and 
U2{r)/[S{S + , respectively. H3 with effective po- 

tential U3{r) is the new interaction arising from random 
variations in the pairing interaction. pl|, p3| In H2, S 
and S are the spin vector and quantum number of the 
magnetic impurities, respectively, and 7? represents the 
quasiparticle spin eigenvector. We assume the spatial av- 
erages of all of the random potentials vanish, {Ui{r)) = 
for j = 1, 2, 3. In the absence of all defects, the order 
parameter Ao(T) = y(c|(r)cj(r)) satisfies the standard 
BCS gap equation, 

Ao = -VT J2 J (2^Tr[p2a2Go(k,c.„)], 

Go"^(k,cj„) = ia;„poo-o - CkPao-Q - Aop20'2, (7) 

where Go is the bare Green's function matrix, y < is 
the uniform (BCS) part of the pairing interaction, iV(0) 
is the single-spin quasiparticle density of states, luq is a 
BCS-like cutoff, and ujn are the Matsubara frequencies. 

We assumed a real bare uniform order parameter Aq, 
and restricted our consideration in H3 to spatial fluctu- 
ations of the amplitude of Aq. The model can also treat 
spatial fluctuations of the phase of Aq by letting AoPicr2 
and U3{y)p2<J2 be generalized to Aoipicr2 + Ao2P20'2 and 
U3i{Y)pia2 + J732(r)P2tT2, respectively 

Our main interest lies in studying H^. Using quan- 
tum Monte Carlo techniques to study a two-dimensional 
square lattice with an on-site attractive Hubbard pairing 
interaction in ffo, [p3| Ghosal et al. obtained interesting 
results in excellent qualitative agreement with those ob- 
tained from STM measurements. p2| We also consider Hi 
and H2 for comparison, because the combination of one 
or both of them with H^ can lead to interesting novel 
behavior. In the Born approximation, these interactions 
add or subtract in a simple fashion. However, in the 
T-matrix approximation for a single defect site, these in- 
teractions interact in a highly non-trivial manner. 

In the self-consistent Born approximation, the quasi- 
particle self-energy E = Si + S2 + S3, where 

S,(k,c^„) = n,^V-(k-k')G(k',c^„)1/,(k'-k),(8) 



Vi (k) , Ui (k) are the spatial Fourier transforms of 
Vi{v),Ui{r), respectively. Neglecting any possible 
anisotropy arising from Fermi surface integrations, the 
effective rates of the three processes are 



1/n = 27rn,7V(0)|C/,(kF)|^ 



(10) 



where rii is the density of defects of type i. 

As in the usual pairbreaking theory, ||l|, ||, |[ |j, G 
has the same form as does Gq, except that ojn and Aq 
are replaced by their renormalized equivalents Cjn and A, 
respectively. We then obtain the standard equations for 
the renormalized gap and Matsubara frequency. 



Cbn = aj„ -I- (1/ri -(- 1/Tpb) 

A = Ao + (1/ri - 1/Tpb) 
1/Tpb = l/r2 + l/r3 



UJr, 



2p2 + A2]l/2' 
A 

2p2 A2]l/2' 



(11) 

(12) 
(13) 



is the total pairbreaking rate. The new physics arises 
from H^. Evidently, within the self-consistent Born ap- 
proximation, the effects of the random interactions are 
exactly equivalent to those of magnetic impurities. 

Using standard pairbreaking theory, H, ||, ||, Q one 
finds 



-— = u(l , 

Ao VTT^ 

U = Wn/A, 

C = l/(Tp6Ao). 

The spatial average gap A(T) is then 

d2k 



(14) 

(15) 
(16) 
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(27r)3 



= 7r|y|7V(0)T 



Tr[p2cr2G(k,a;„)], 
1 



7^ vT 

aj„ I <Wo 



leading to the standard equation for Tc/Tcq ~ t, 

■1 



= ln(t) 
ttpb = l/(T-pbT'co), 



^^2 2TTtJ 



(17) 



(18) 
(19) 



G-i(k,c^„) = G(7i(k,c^„)-S(k,c^„), 



(9) 



where ijj{x) is the digamma function. For small apb, 
Tc ~ Tco — 7r/4rpb. We note that Tc/Tco can be sup- 
pressed to zero even without any magnetic impurities, 
for l/r3 > l/r3c — where 7 = 1.781 is the ex- 

ponential of Euler's constant. In addition, the supercon- 
ductivity becomes gapless for 1 > t^c/t^ > \2e~'"l^ ~ 
0.912. Thus, even an isotropic, s-wave superconductor 
can become gapless, as observed in the cuprates with 
STM.|,[lO[ ^ 

This can only occur in short coherence length super- 
conductors with strong local inhomogeneities in the pair- 
ing interaction, as is a likely explanation for the vanishing 
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of the Tc in the highly underdoped region of the cuprate 
phase diagram, although that region is also complicated 
by the simultaneous appearance of local SDW/CDW or- 
der at the non-superconducting regions not included in 
this calculation, 

In order to make direct comparison with STM experi- 
ments, we solve the T-matrix for a single defect site. We 
assume the site has all three types of defects associated 
with it. We approximate the effects of the magnetic im- 
purity by assuming that its spin behaves classically. [p4| 
Then, the T-matrix equation can be solved exactly, 



V{0) 



f(uJn) = -. 

l-go{iOn)V{0) 

vio) = Y^vm, 

i=l 

-ttN{0) 



(20) 
(21) 



So(w„) 



A2]i/2 



[iw„po(To -f A0P2CT2], (22) 



is the bare Green's function at the origin, 1 ~ Pqctq is the 
rank 4 identity matrix, and the effects of a finite quasi- 
particle energy bandwidth have been neglected. Bound 
states within the gap at T = at the frequency lu are 
obtained from 



det\i-goiiu;)Vm = 0. 



(23) 



Solving Eqs. ( pT| ) - ( p3| ) exactly, we generally obtain 
four bound states within the gap at a; = cUf,Ao, where 



LJb = 


±[A + sRYI^/B, 


(24) 


A = 


IQvlvl ^ a\{al_ -l-4f^). 


(25) 


B = 


al + Avl, 


(26) 


C = 




(27) 


R = 




(28) 


a± = 


v\ + vl-vl± 1, 


(29) 


Vi = 


7r7V(0)?7,(0), 


(30) 



and s = ±1. 

It is useful to rewrite Eq. ( p^ for the three special 
cases of two defects only. For non-magnetic defects, V2 = 
0, there are only two bound states symmetric about zero 
bias, 



± 



1/2 



(1 



(31) 



For the trivial case of a site with just a non-magnetic 
impurity, Eq. (|3l| ) shows that there are no poles inside 
the gap. However, for a pairing interaction defect alone, 
there are poles at ZUb = ±|1 — -I- f|). 

For impurities alone, W3 = 0, there are also only two 
bound states symmetric about zero bias, at 



± 



|1 -|- U+U- 



[(1-f 0(l+u2_)]l/2 



(32) 
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FIG. 1: Sketch of a bound state at zero bias with v = 
(0,0, ±1), along with data from Ref. (9) for the STM dif- 
ferential conductance above a Zn site in BSCCO. 



where u± — ±U2. For v\ — 0, this reduces to the result 
of Shiba for a classical magnetic impurity, cUb = ±|1 — 



even though we did not average over the 



classical spin direction before summing the T-matrix. |24| 
The most interesting case arises when both U2 , W3 7^ 
and U3 and ^2 7^ i^3- In the limit wi = 0, there are four 
bound states symmetric about zero bias at 



± 



{v. 



(H-t;2)(l + t,2) 



(33) 



where s = ±1, f -t = V2 ± fs- If U2,W3 7^ and vi 7^ 113, 
there arc four bound states, exhibiting reflection symme- 
try about zero bias. Otherwise, if either vi or W3 — 0, 
or if ^2 = ^3, there are only two bound states that are 
symmetric about zero bias. We note that Eq. (|33| ) for 
either t;2 = or U3 = reduces to Eq. (^Tj) and (|32|) with 
v\ = 0, respectively. In addition, for V2 = ^3, it reduces 
to Eq. (IH) with v\ — ^ and vi — > 2w2, etc. More- 
over, setting V2 ~ W3 in Eq. ( p^ ) with v\ arbitrary leads 
to only two bound states symmetric about zero bias, at 
W6 = ±[|5_|/6+]i/2, where 6± = {\+vlY±Avl. Thus, we 
conclude that in order to obtain four bound states, two 
on each side of zero bias, one requires 7^ W2 7^ ?^3 7^ 0. 

When the defect is a quantum spin with a single com- 
ponent normal to the surface, the spin operator com- 
mutes with the Hamiltonian, and the spin states are 
easily described by \SM), with S^\SM) = M\SM). 
Then the magnetic impurity in the presence of the non- 
magnetic potential and the pairing disorder can all be 
solved exactly. There are bound states for each of the 
25-1-1 eigenstates. 

In Figs. 1 and 2, we have illustrated how this so- 
lution can aid in understanding the STM results ob- 
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site Hence, a generalization of our results to d- 

wave superconductors would not be straightforward. 

In summary, we have shown that disordered short co- 
herence length superconductors can exhibit pairbreaking 
from spatial fluctuations in the pairing interaction, in a 
manner very similar to that found with magnetic impu- 
rities. We studied the effects of a single site with up to 
three types of defects using the T-matrix approximation, 
and found bound states within the superconducting gap 
arising from either pairing fluctuations or magnetic impu- 
rities. Our best fits to the scanning tunneling microscopy 
data above the sites of Zn and Ni impurities suggest that 
Zn behaves either as a strong pairing fluctuation defect 
or as a strong magnetic impurity, whereas Ni behaves as 
both a strong magnetic impurity and a strong pairing 
defect. 



FIG. 2: Sketch of the two positive bias bounds states obtained 
with V = (0, ±0.51, ±1.01), along with the data of Ref. (10) 
obtained above the Ni site in BSCCO. 



tained from the Bi sites directly above Zn and Ni im- 
purity sites in the presumed top underlying Cu02 plane 
of BSCCO, respectively.!, |l^ In Fig. 1, we fit the Zn 
STM data with a single pole obtained from Eq. ( p^ ) 
with V = {vi,V2,V3) = (0,0, ±1) [or equivalently with 
V ~ (0,±1,0)] and a width Sco chosen to fit the data. 
The experimental peak center appears at a slight offset 
from zero, which can be understood quantitatively by 
adjusting wi << 1 and a width that merges the two re- 
sulting peaks into a single peak broader than the offest. 

In Fig. 2, we used Eq. @) with v = (0, ±0.51, ±1.01) 
[or equivalently with v — (0, ±1.01, ±0.51)] to fit the two 
peak positions, and adjusted the widths to those of the 
data obtained above a Ni site.[|l^] We did not show the 
other two peaks that one expects from Eq. (^3|). That 
is because the data not pictured here show that the two 
peaks present in these data appear at equivalent negative 
biases on adjacent Bi sites, [0 for reasons that are not 
particularly clear. In any event, our theory suggests that 
the STM data for Zn are consistent with it behaving ei- 
ther as a strong pairing fluctuation defect or as a strong 
magnetic impurity, and the Ni data suggests it behaves 
as both a strong magnetic impurity and a strong pair- 
ing defect in BSCCO. The similar strengths of the two 
pairbreaking defects might explain why the Tc suppres- 
sions obtained in doping BSCCO with these elements are 
nearly identical. 

Thus, we solved the T-matrix in this modified BCS 
model of a local, on-site attractive pairing interaction 
with three types of defects on a site.Q For a supercon- 
ductor with local, near-neighbor pairing of d^2_j^2-wave 
symmetry, the local gap A^- at the site («, j) on a tetrag- 
onal lattice is coupled heirarchically to the Ai^ji at every 
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